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fraction  to  the  ratio  of  two  contiguous  Bessel  functions «  is 
exploited  to  derive  the  representation  of  a  single  Bessel 
function  in  terms  of  a  series  of  products  of  continued  fractions. 
This  representation  is  then  made  the  basis  of  a  new  method  for 
computing  niimerloally  Bessel  functions  of  the  first  kind.  The 
method  is  seen  to  be  especially  valuable  in  computing  ^ssel 
functions  of  order  much  larger  than  unity. 


PUBLISHED  DECEMBER  1963 


U.  S.  NAVAL  ORDNANCE  LABORATORY 
WHITB  OAK,  MARYLAND 


i 


NOLTR  63-228 


16  October  I963 


A  NOTE  ON  BESSEL  FUNCTIONS  AND  CONTINUED  FRACTIONS 


The  work  reported  herein  was  carried  out  In  RB  under  Foundational 
Research,  FR-43. 

The  contents  of  this  report  will  be  Incorporated  In  an  article 
concerning  hyperge<»ietrlc  functions  and  continued  fractions  which 
Is  to  be  submitted  for  publication  In  the  open  mathematical 
literature. 

This  document  Is  Intended  to  serve  as  a  progress  report. 

The  author  Is  grateful  to  Ilr.  A.  H.  Van  Tuyl  for  several 
discussions,  and  to  Susan  Nadlgosky  for  programming  the  machine 
calculations . 


R.  E.  OEBNINO 
Captain,  USN 
Commander 

Z.  I.  SLAVSKY 
By  direction 


11 


NOLTR  63-228 


CONTENTS 

Pag* 

INTRODUCTION . 1 

BESSEL  FUNCTIONS  AND  CONTINUED  FRACTIONS .  1 

CONVERGENCE  QUESTIONS .  4 

EXAMPLES  OF  NUMERICAL  COMPUTATION .  7 

ILLUSTRATIONS 

Figure  Title  Page 

1  Values  of  M  such  that  J2m(^)  ^  .  ^ 


NOLTR  63-228 


INTRODDCTION 


1.  The  purpose  of  this  report  is  twofold.  First  we  will 
derive  several  examples  of  a  class  of  not  widely  known  Identities 
which  relate  Bessel  functions  to  certain  contAaued  fractions. 
Second,  with  the  aid  of  one  of  these  Identities  we  will  develop 

a  new  method  for  the  numerical  generation  of  Bessel  functions  by 
hl^  speed  digital  computers.  Throughout  the  report  our  Interest 
will  center  on  Bessel  functions  of  the  first  kind,  Jn(z)«  ^are 
n  Is  any  Integer  and  z  Is  an  arbitrary  real  number,  although  It 
will  be  apparent  that  generalizations  to  other  cylinder  func¬ 
tions  of  Integer  order  and  arbitrary  complex  argument  are 
possible. 

2.  As  for  the  nximerlcal  computation  of  Bessel  functions, 

we  have  In  mind  not  the  preparation  of  tables  but  the  calculation 
of  Jn(a)  In  the  midst  of  a  lengthy  sequence  of  calculations  by  a 
digital  computer  where  the  value  of  z  Is  not  known  ahead  of  time. 
The  greatest  need  for  such  a  routine  arises  when  z  and  n  are  both 
large.  When  n  »  z,  the  computation  of  Jn(z)  may  be  efficiently 
carried  out  with  the  series  representation,  and  when  z  »  n,  the 
asymptotic  expansion  provides  the  most  direct  method.  However 
there  Is  no  completely  reliable  standaz^  method  (to  our  knowledge) 
for  computing  Bessel  functions  In  the  Intermediate  region  where 
z  ^  n,  especially  when  n  >  20.  To  check  the  accuracy  of  the 
method  developed  here  several  short  tables  of  Bessel  functions 
have  been  made  and  Included  In  this  report.  In  no  ease  has  the 
6th  decimal  place  been  discovered  to  be  Incorrect. 

BBSSSL  FUNCTIONS  AND  C(»ITINUBD  FRACTIONS 

- 3.  Since  the  time  of  Bessel  hlmseTr  It  has  been  known  that 

the  ratio  of  two  contiguous  Bessel  functions  may  be  written  as 
an  Infinite  continued  fraction.  *nils  fact  follows  In  an  ali¬ 
mentary  way  from  the  recurrende relations  (Watscm,  Ibid.,  pp.  46, 
133 K  which  may  be  rearranged  In  the  form: 

Jn(»)  ^  I  H  (1) 

Jn-l(*)  *  ,  Jn->-l(^ 

*2  nJ„(z) 


•A  cylinder  function  Is  any  function  satisfying  the  Bessel 
fimctlbn  recurrence  relations.  For  a  full  description,  see 
Watson,  Theory  of  Bessel  Functions  (Caidbrldge  Unlv.  Press,  193o)» 

p.  82. 
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« 


The  extension  of  this  Identity  into  explicit  continued  motion 
form  is  obvious.  If  we  chsnge  traia  Jn(z)  to  Jn(2x)  to  avoid 
factors  of  1/2  throughout «  we  obtain 


Jn(2x)  ^  I _ 

Jn-l(2x)  ^  _  xVn(nfl) 

1  -  . 


(2) 


In  the  shorthand  notation  which  we  will  \ise  throughout,  Bq.  (  1  ) 
is  simply 


Jp(^)  .  _5l. 

J„.l(2x)  1-Pn 


(3) 


where  Xn  •  ^ ,  and  Fh  ie  the  specific  continued  fraction  defined 
for  all  n  1  by  the  functional  relation 


«  Xnxn4-1  (4) 

Clearly  then,  it  is  possible  to  represent  the  ratio  of  amr  two 
Bessel  functions,  J»fk(2x)/Jai(2x),  idiose  orders  differ  iy  an 
integer,  by  a  product  of  continued  fractions t 


Such  a  product  representation  of  the  ratio  of  two  Bessel 
functions  we  will  denote  Fi'^K(2x).  We  also  define  P|  9  1;  and, 
if  n  >  P®  »  l/PJ. 

4.  With  these  few  preliminaries  taken  care  of,  we  now 
state  the  identity  which  we  will  use  in  the  numerical  oemnuta- 
tionsf  paiagnacftt.lO  •  Several  generalisations  of  this  identity 
will  be  given  in  the  following  paragbaph.  the  identity  is; 
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1 

Jn(*) 


W 


>2bi 

N 


(*) 


maO 


where  €©  *  1*  €  m  •  2  for  m  >  |  . 

A  proof  of  Eq.  (  6  )  Is  arranged  very  easily.  It  depends 
trivially  on  the  fact  that  the  P's  are  ratios  of  Bessel  functions^ 
and  the  well-known  addition  theorem  (Watson,  Ibid.,  p.  23) t 


J2m(*)  »  1*  (7) 

•-0 

Thus  we  have  arrived  at  '  new,  or  at  least  not  widely  known, 
result: (*)  the  representation  of  a  single  Bessel  function, 
rather  than  the  ratio  of  two  Bessel  functions,  by  Infinite 
continued  fractions. 


5.  It  should  be  mentioned  that  other  analogous  continued- 
fraction  representations  of  Bessel  functions  can  be  easily 
arrived  at.  For  example.  In  quantum  elect rodyammlcs  the  sunma- 
tlon  of  certain  classes  of  Feynman  diagrams  can  be  shown (**)  to 
lead  to  the  following  Identity: 


In  idilch  the  square  of  a  Bessel  function  Is  given  by  a  series 
of  squares  of  products  M  (z).  This  representation  depends  on 
the  well-known  Identity  (Watson,  Ibid.,  p.  338): 

S  .1.  (9) 

Also  we  have: 


1  a  a  S  (-1)"  .  (10) 

Jn(s)  sin  z  m^ 

(•)to  the  author's  knowledge  the  relations  given  In  Bqa.  (6^8.,)1<>^)  <iD 
not  appear  In  the  published  literature  an  Bessel  functions.  ^ 
view  of  the  simplicity  of  derivation,  however.  It  seems  unlikely 
that  they  have  heretofore  gone  unnoticed. 

f**)Z.  Fried  and  J.  H.  Sberly  (paper  In  preparation).  It  was 
this  work  lAiloh  led  the  author's  attention  to  the  subject  of  the 
present  report. 
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and  several  others  slightly  more  oomplloated^  idiloh  we  onit  here* 
but  each  of  wh&oh  Is  related  in  an  obvious  way  to  a  known  tuama- 
tlon  Identity  like  Eq.  (  7  )  or  (  S  ). 

6.  Some  similar  results  appear  to  be  possible  with  the 
Bessel  funetlons  replaced  by  general  confluent  hypergeosMtrlo 
functions.  However  these  generalizations  are  peripheral  to  our^ 
purpose  here,  and  we  will  restrict  the  present  dlseusslon  to 
Bessel  functions  of  Integer  order. 

CONVroOENCE  QUESTIOMS 

Tii  A  very  significant  property,  fraa  the  standpoint  of 
numerical  computation,  of  the  sums  like  Eq.  (  6  )  le  their  quite 
rapid  convergence.  Before  going  further  Into  the  suitter,  though, 
several  points  should  be  made  clear.  In  the  first  place,  all  of 
the  Infinite  series  and  Infinite  continued  fractions  which  will 
be  dealt  with  do  converge  (Watson,  Ibid.,  pp.  133-134,  et  passim). 
The  series  and  the  continued  fractions  In  all  oases  represent 
analytic  functions  of  both  Index  and  argument,  with  simple  poles, 
at  worst,  for  singularities.  That  Is,  the  question  of  conver¬ 
gence  Is  settled  from  the  outset,  and  our  concern  Is  entirely 
for  the  practical  question  of  speed  of  convergence.  Specifically, 
we  will  be  Interested  In  the  answers  to  these  two  mcstl<»sz 
Given  the  continued  fraction  representation  for  J||(s),  Eq.  (  6  )• 
how  many  terms  of  the  Infinite  series  must  we  keep  to  ensure  the 
desired  nxuaerlcal  accuracy;  and  how  shall  we  truncate  the 
Infinite  continued  fractions  appearing  In  the  terms  of  the 
series?  Thus  we  are  Intetrested  In  making  sufficiently  accurate 
approximations  to  the  terms  In  the  series  as  well  as  to  the 
series  Itself.  We  shall  answer  the  first  question  first.  Assume 
that  we  want  to  compute  Jn(z),  end  that  n  and  z  are  both  large 
and  comparable  im  magnitude  so  that  neither  the  usual  series 
representation  nor  the  asymptotic  expression  for  ere  conven¬ 
iently  usable.  Thus  we  turn  to  the  series  (  6  ).  From  the 
definition  of  the  P's  in  terms  of  ratios  of  Bessel  functions  It 
Is  clear,  though,  that  as  far  as  speed  of  convergenoe  goes  we 
need  consider  only  the  problem  of  making  N  large  enough  so  that 
S||  Is  sufficiently  close  to  unity,  where  Sn  is  defined  to  be: 

M 

Sm(z)  ■  ^fenJ2n(*) 

We  can  conclude  that  convergenoe  Is  likely  to  be  quite  rapid 
after  we  reach  the  terms  2n  ^  z  from  Inspection  of  the  inequality 
(Watson,  Ibid.,  p.  49): 
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ft,  in  mmabAnm  sunmatlon.  In  which  the  computer  Is  usually 
limited  to  eight  significant  figures >  it  is  of  course  useless 
to  consider  terms  which  are  smaller  than  lO'O  times  the  largest 
term  in  the  series «  so  we  are  Interested  in  knowing,  for  a 
given  z,  how  big  the  J2m(z)  get.  Although  the  maximum  value  of 
J2m(z)  is  known  to  occur  near  the  point  z  »  2m,  the  problem  of 
determining  boxuids  on  the  size  of  the  maximum  seems  to  be 
unexplored  territory.  Therefore  we  will  be  satisfied  with  the 
crude,  but  almost  certainly  true(*),  estimate  that  for  any 
z  <  100  there  will  be  a  term  in  the  series  Ses,  Eq.  (  7  ), 
greater  than  0.10.  This  means  that  we  need  take,  in  the  series 
Sn,  M  only  so  high  that  J2p(z)  <  10" 9  for  all  p  >  M.  A  glance 
at  tables  of  Bessel  functions  (**)  enables  us  to  construct  a 
graph  to  serve  as  a  guide  in  choosing  this  value  of  M,  given  a 
value  of  z.  Roughly  speaking,  the  graph  shows  that  a  sufficient 
value  of  M  is  given  by  M  «  ^  z  Hr  6.  The  graph  is  reproduced 
below  in  Figure  1. 

Now  we  must  attend  to  the  second  question,  that  of 
approximating  the  infinite  continued  fractions  which  appear  in 
every  term  of  the  series  (  6  ).  Notice  that  by  repeated  use  of 
the  functional  relation  (  4  )  we  may  easily  compute  all  of  the 
Fn  if  we  already  know  one  of  them.  Also  notice  that  the  higher 
the  subscript  n,  the  more  rapidly  F^  converges.  Thus  our  pro¬ 
cedure  is  vex*y  simple:  we  determine  what  the  highest  occurring 
subscript  is,  compute  the  corresponding  F  from  its  form  as  a 
continued  fraction,  and  then  compute  all  other  F's  from  it  by 
means  of  the  relation  (  4  ).  The  computation  of  the  highest  F, 
say  Fp  ,  is  straightforward.  Written  out  explicitly,  Fp  is: 


Fp  •  _ _ 

1  .  Xp+lfEiS 

1-  ^PH-3XpH-4 
1-  . 


The  first  step  is  to  comoute  only  the  first  level,  XpXp^j* 
is  already  less  than  10"”,  then  F  •  0  is  a  sufficiently  accurate 
approximation.  If  not,  then  the  second  level  is  included  and  the 
fraction  computed  to  this  approximation.  If  the.inolusion  of  the 
second  level  makes  a  difference  of  less  than  10"^  in  the 

(*)This  estimate  was  arrived  at  by  inspection  of  the  Tables  of 
Bessel  Functions  of  the  First  Kind  (Harvard  Ihiiv.  Press,  1947- 
1951).  Here  are  tabulated  Jq  to  J133  for  0  <  z  <  100. 

(**)Very  convenient  for  this  purpose  is  a  table  in  Jahnke  b  Bmde, 
Tables  of  Function^,  4th  ed.  (Dover  Publ.,  NCw  York,  1945)# 

pp.  171-179. 
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single-level  result,  then  Pp  ■  XpXp+i  is  good  enough;  and  so  on 
until  the  inclusion  of  one  more  level  affects  the  preceding 
result  by  less  than  10-°,  The  only  remaining  matter,  the 
determination  of  the  highest  occurring  subscript,  is  easily 
settled.  If  we  are  computing  Jm(z),  and  have  decided  to  keep  M 
terms  in  the  series  (6  ),  then  the  larger  of  N  and  2M  will  be 
the  largest  subscript. 

EXAMPLES  OF  NUMERICAL  COMPUTATION 

10.  On  the  following  pages  we  display  some  results  obtained 
by  machine,  computation  of  Bessel  functions  using  the  continued 
fraction  method  described  in  this  report.  The  labeling  should 
be  clear.  For  a  given  value  of  z  several  different  Jn's  are 
computed;  then  z  is  changed  and  the  J)|'s  are  computed  again, 
and  so  on.  Most  of  the  values  have  been  checked  so  far  as 
possible  against  standard  tables,  in  particular  the  Harvard 
tables  mentioned  in  Paragraph  8  ,  and  no  errors  have  been  dis¬ 
covered  in  the  sixth  decimal  place.  Only  near  the  zeros  of  the 
J])}(z)  are  there  errors  in  the  seventh  decimal  place.  In  at 
least  one  case,  Jo5(10),  there  was  no  error  until  the  twenty- 
fourth  decimal  place. 
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Z  =  30 


J36 

0.7l29A3aflE-C2 


J24 

-0.32381186E-CI 


J12 

C. 


Z  =  29 


J36 

0.3437l736e-0? 


0.789852976-01 


J12 

0. 1199R265E-CC 


Z  ^  28 


436 

0.15<i76065e-C2 


J26 

C.17043918E-0C 


J12 

-0.38292769E-C.' 


Z  -  27 


J36 

0.669f6067E-03 


J26 

0.22Ca7C99E-CC 


J12 

-D.12946655E-00 


J3C 

C. 16393584  ;>C0 


Ji8 

0.15647528L-00 


JZ 

0.48622736€-02 


J3C 

0.lC304804r-00 


J18 

0.141372l3f-C0 


46 

U. 12702166L-C0 


J20 

0.676853836-01 


J18 

0.39404874-01 


J6 

0.13933615£-00 


J3'' 

0.409592266-01 


418 

0.89281379'.-01 


46 

0.271277286-01 
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L  *  26 


J36 

0.25380371E-C3 


J24 

O.227l43«7E-0C 


J12 

-0.161090376-00 


2  »  25 


J36 

0.91988669E-OA 


J2A 

0.1997785CE-C0 


J12 

-0.72867787E-01 


2  -  2A 


J36 

0*30830362E-04 


J2A 

O.155OA220E-CC 


J12 

0.7299011AB-01 


2-73 


J36 

0.99162S13E-05 


J2A 

0. 107822386-00 


J12 

0.173018576-00 


J30 

0.22882923F-01 


J18 

-C.17521880F-00 


J6 

-0. 11366042F-00 


J10 

0.ll809026!:-0l 


J18 

0.17377267L-C0 


J6 

-0.15870029F-00 


J30 

C.56256811L-02 


J18 

-C.93111813C-01 


J6 

0.6A5A7005F-01 


J30 

0,?A69772ir-02 


JIS 

C.3AOl90A7t-0l 


J6 

0.90859228F-01 
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1  ^  ?:i 

J16 

0.2692  1725E-C‘J 

J30 

0.99659806--C3 

J2« 

C.67729434E-01 

JIP 

U.15492100I--C0 

J12 

0.l5657e7lE-C0 

Jt 

0.  l7325246t:-00 

Z  »  21 

ja6 

0.69409896E-C6 

J30 

C.36e22633}>03 

J2A 

C.38570C37E-01 

J18 

0.23l!>726lf-00 

J12 

a.32928689E>0l 

J6 

C. 1076486lh-00 

» 

L  2  - 

J36 

C.16:00121E-C6 

J3C 

0.124».l536r-03 

J29 

0.1992910SE-01 

J18 

0. 25108984. -CO 

J12 

-0.11899C69E-CC 

J6 

-0.55.fl6067.--0l 

Z  -  19 

1 

J36 

0.33960712E-C7 

J2C 

0.37849l46r-04 

\ 

\ 

J24 

C.93306652E-C? 

J18 

0.223523l4f>C0 

\ 

J12 

-0.20545822E-CC 

J6 

-0. l78767lSt-00 

10 

NOLTR  63-228 

i  »  18 

J16 

0.63352702E-C8 

J30 

C.10393653C-04 

0.3945fll32E-C2 

J18 

0.17062988L-00 

Jt2 

-C.17624115E-00 

J6 

-0.155956191-00 

1 

Z  »  17 

J36 

0.10601692E-C8 

J30 

C.25460C67F-C5 

J2M 

C.16996625E-C2 

Jl« 

C.ll381l0ir-00 

J12 

“O.60574816E-OI 

J6 

0.7 1536081 F-03 

/  «  16 

J36 

C.14e35t78E-C9 

J3C 

0.55052391(^-06 

J24 

0.50874504E-CJ 

Jl« 

C.6684808lf;-Cl 

1 

1 

! 

i 

J12 

0.11240C26E-OC 

J6 

0.16672073^-00 

I 

! 

i 

1 

Z  »  15 

J36 

0.18091641E-IC 

J30 

0.103747liL-C6 

i 

\ 

1 

J26 

0.19266956E-03 

J16 

0.346259820-01 

1  < 

r 

1 

1 

i 

J12 

0.23666584E-CC 

J6 

0.20614970L-00 
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z  » 


J36 

0.l85072l8e-ll 


J24 

0.400638926-04 


J12 

C.28545026E-0C 


Z  »  13 


J36 

C. 155120686-12 


J24 

0. 906046346-05 


412 

0.261536876-00 


Z  »  12 


J36 

0.103455786-13 


424 

0.173322636-05 


412 

C.19528018E-CC 


Z  »  II 


436 

0. 529085686-15 


424 

0.273628116-06 


412 

0.12159979E-C0 


J30 

C.16775401"-0T 
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